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Abstract 

The process of cascaded downconversion and sum-frequency generation inside an optical cavity 
has been predicted to be a potential source of three-mode continuous- variable entanglement. When 
the cavity is pumped by two fields, the threshold properties have been analysed, showing that these 
are more complicated than in well-known processes such as optical parametric oscillation. When 
there is only a single pumping field, the entanglement properties have been calculated using a 
linearised fluctuation analysis, but without any consideration of the threshold properties or critical 
operating points of the system. In this work we extend this analysis to demonstrate that the singly 
pumped system demonstrates a rich range of threshold behaviour when quantisation of the pump 
field is taken into account and that asymmetric polychromatic entanglement is available over a 
wide range of operational parameters. 

PACS numbers: 42.50.Dv,42.65.Lm,03.65.Ud,03.67.Mn 
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I. INTRODUCTION 



The modern field of quantum information originally focused on what is known as 
discrete- variable entanglement and developed to include the study of entanglement between 
continuous-variable phase quadratures of the electromagnetic field, which have a close anal- 
ogy with the original position and momentum considered in the famous Einstein-Podolsky- 
Rosen (EPR) paradox [1] . In the beginning this research considered bipartite entanglement 
as produced by, for example, the optical parametric oscillator (OPO) [2J and led to experi- 
mental demonstrations of the EPR paradox [3| and of what is known as continuous- variable 
quantum teleportation . Many systems have now been studied, both theoretically 

and experimentally, with continuous-variable bipartite entanglement now considered an im- 
portant resource for quantum information applications 7J. 

Recently there has been much attention paid to the production of continuous-variable 
tripartite entanglement, obtained either by mixing squeezed beams on unbalanced beam- 
splitters js], 9], or via the interaction of multiple input beams in nonlinear media with cas- 
or concurrent y£j* nonlinearities. Among the latter are systems using either single 10] , 



cadec. 
twin 
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current Y_> 



or triple 
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18] nonlinearities. These nonlinear processes have 



been analysed and demonstrated in both the travelling-wave and intracavity configurations. 
In this work we are interested in an intracavity process which combines parametric down- 
conversion with sum frequency generation, as theoretically analysed by Yu et al. 19j. The 
idea of combining these two processes is due to Smithers and Lu [201], who did not consider 



enclosing the processes in an optica 
was first analysed by Guo et al 



cavity. The intracavity process with two pump fields 
121 ]. who used quantum Langevin equations 21| with an 
undepleted pump approximation, which can give no insight into threshold behaviour or any 
critical operating points. An analysis which included quantisation of the two pump fields 
was performed by Olsen and Bradley [ljj, demonstrating that the system had quite dif- 
ferent stability and threshold behaviours to the normal OPO. As Yu et al. [l9] have also 
used an undepleted pump approximation (also known as the parametric approximation) and 
quantum Langevin equations for the singly pumped intracavity system, they are not able 
to determine the threshold behaviours or the stability of the equations they use. In this 
paper we apply a fully quantised treatment of all the interacting fields, finding that there 
are two separate parameter regimes, one of which has an oscillation threshold while the 
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other does not. The reason for this behaviour, which is more complicated than that of the 
standard OPO, is that downconversion considered separately does exhibit a threshold while 
sum frequency conversion does not. As we will demonstrate below, it is the competition and 
interplay of these two processes which leads to more complicated behaviour. 



II. HAMILTON! AN AND EQUATIONS OF MOTION 

Our system is a nonlinear medium inside an optical cavity which is pumped at frequency 
ujo and is resonant at all the frequencies involved. In Yu et al. [l9( the nonlinear medium is a 
quasiperiodic superlattice. In the downconversion part of the intracavity process, two fields 
at u)\ and uj% are generated, where uj = U\ + uj 3 . We will denote the effective nonlinearity 
for this process by xi- The pump field at ujq can then combine with the field at 003 in a sum 
frequency generation process, to produce a further field at u>2, with the effective nonlinearity 
represented as \2- We will use the annihilation operator b to describe the field at ujq, while 
the operators dj will be used for the fields at ujj. 

The Hamiltonian can be written as 

Htot — H pump + H\ nt + Hd&mp, (1) 

where the interaction Hamiltonian is 

-ffint = ihxi{ba\al - 6 f aia 3 ) + ihx 2 {ba 3 al - tfa\a 2 ), (2) 
the Hamiltonian describing the cavity pumping is 

H m = th(etf -e*b), (3) 
and the cavity damping Hamiltonian is 

= HT 6 f + f $S) + h £ (fj-at + fta,-). (4) 

j=l:3 

In the above, e is the pump field which enters the cavity, which will be described classically, 
and the Tj are reservoir operators for each of the intracavity modes. 

To calculate the fluctuation and entanglement properties of the system we will derive 



fully quantum equations of motion using the positive-P pseudprobability distribution [221 ]. 



as this naturally allows us to calculate the normally-ordered operator expectation values 



required to find output spectra. Proceeding via the normal methods 23j and making the 



zero temperature and Markov approximations for the reservoir 
Planck equation for the P function 25|, |26) of the system as 
dP { 



24j . we find the Fokker- 
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P, (5) 



K daida 3 da{da* 3 J \da 5 d/3 da* 3 d/3* 

where the jj are the cavity loss rates at frequency ujj. As this Fokker-Planck equation does 
not possess a positive-definite diffusion matrix, we must double the phase-space and use 
the positive-P representation to find the appropriate stochastic differential equations. This 
results in replacement of the conjugate variables by and (3 + , which are the complex 
conjugates of the uncrossed variables only in the mean. Stochastic averages of products of 
these variables are then equal to normally-ordered expectation values of the corresponding 
operators. We find the coupled set of stochastic differential equations as 
da.\ 



'XiP 
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where the rjj are real Gaussian noise terms with the properties 



7fe(t) = 0, r ] ]{t)Vk{t')=5 jk 5{t-t'). (7) 

In cases where the procedure is valid, the noise terms may be dropped and the resulting 
semiclassical equations linearised about their steady states, which results in the process be- 
ing treated as an Ornstein-Uhlenbeck process [27], allowing for easy calculation of the output 
spectra. The validity of this linearised fluctuation analysis is usually found by calculating 
the eigenvalues of the resulting drift matrix for the fluctuations and requires knowledge of 
the classical steady-state solutions. In fact, in the present case, we find that stochastic 
integration of the above equations presents various stability problems in the regions where 



they cannot 



sentation 



De linearised, so that in section [VI] we will turn to the truncated Wigner repre- 



28] to find time domain solutions in these parameter regimes. 



III. LINEARISED FLUCTUATION ANALYSIS 

In the steady-state, we can always decompose the system variables into their mean values 
and a part which fluctuates about these. In many cases the mean value solutions of the 
noiseless equations are equal to the operator expectation values and the fluctuations can 
be treated as being stable and Gaussian about zero means. In these cases we may use 
this linearised fluctuation analysis as a simple method to calculate measurable spectra. We 
will perform this process on the positive-P equations, beginning with the decomposition 
on = Wi + 5ai and similarly for (3. This gives us the set of equations for the fluctuating terms 



San = -7i5ai + XiP^at, + \ ^-(vi +^2), 



dt 



dt 

— 5a* 2 = -<y 2 5a* 2 + X2p*5a* 3 , 



d 5a 3 = -j 3 5a 3 + Xx^a* x - X2p*5a 2 + \ ^-{Vi - ^2), 



dt 



j t 6/3 = -7tM 
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j t 5(3* = - lo 5[3*, (8) 
which may be written in matrix form for the vector 

5a = [5ot\, 5a{, 5a 2 , 5a 2 , 5a%, 5al, 5(3, 5(3*] T , (9) 

as 

^-§a — —A 5a + B dW, (10) 
at 

where A is the drift matrix, B contains the steady-state coefficients of the noise terms, and 
dW is a vector of Wiener increments. The condition for stability of the fluctuations is that 
the eigenvalues of A have no negative real parts. When this condition is fulfilled, we may 
calculate the intracavity spectral matrix as 

S{uj) = {A + iul)' 1 BB T (A T - iwtY 1 , (11) 



which, along with the well-known input-output relations 2l| . allow us to calculate the 
measurable spectral quantities outside the cavity. 



A. Steady-state classical solutions 

The classical equations for the mean values are found as 

da~i 



dt 
da 2 

~dT 

da 3 

1t_ 
dp 
~dt 



-7i«i + Xi&W, 
-72 a 2 + X2a 3 P, 

-73«3 + Xl®lP ~ X2«2/3*, 

e - 7o/3 - Xittia3 - X2a 2 a* 3 , (12) 



and may be solved for the steady-state solutions which enable us to perform the necessary 
stability analysis. We find that the solutions are divided into two different classes, depending 
on whether an oscillation threshold is present or not. 
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1. Regime with threshold 



If Xi72 > X2715 we fi n< i that the system has a threshold pumping value below which it 
will not oscillate. If the value of the pump field e is below 



(13) 



71 72 

the signal modes will not be macroscopically occupied. We note that this is totally different 
from the expression which would be expected if we considered the threshold for the down- 
conversion process in isolation, which would be given as e° = 7n^/7i7.s/ Yi ■ This difference 



in threshold cannot be calculated in the approach taken by Yu et al. 19j. The analytical 



expressions for the different mean-value steady-state solutions are (note we will now drop 
the bar over the variables for notational convenience) : 
(i) e < e c 
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(15) 



where 6 is an undetermined phase. We notice that due to the presence of the square root, 
the sign of these solutions is unknown. However, inspection shows that the square roots all 
have to be the same sign, whether this is positive or negative. Nevertheless, because the 
only phase we know is the phase of the pump field e which we take as real, (3 SS will also 
be real. As 9 is not fixed, these above threshold solutions will exhibit phase diffusion, as 
previously found in the nondegenerate parametric oscillator 



2. Regime without threshold 

We find that if X1I2 < x\l\ there is no threshold predicted by the classical equations, 
which means that for any value of the pump field the signal modes will not be macroscopically 
occupied. The expressions for the steady-state solutions are the same as the below threshold 
solutions of the previous case, 

a ss = -, 

7o 

a]' = 0. (16) 

In our analyses in this regime we will scale the pump amplitude by the normal OPO thresh- 
old, e° = 70^/7173/xi, so that the pump field e will be expressed as a proportion of this 
threshold. Now that the classical steady-state values in the different areas are known, we 
can analyse the stability of the fluctuations. 



B. Stability analysis 

To determine the validity of the linearisation process we will now analyse the eigenvalues 
of the drift matrix A of Eq. [TUJ This is written out in full as 
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(17) 

The analytical expressions for the eigenvalues of the matrix A above threshold are not 
easily obtained, but below threshold where all the a s j s are zero and (3 SS = e/70, we find the 
characteristic polynomial as 

(7o - A) 2 [(71 - A)( 72 - A)( 73 - A) + \r\xl ~ Xl) + P ss \lixl ~ 72X?)f = 0, (18) 



S 



By studying the variation of this function we note that the system is always stable below 
threshold (whether X1I2 > X2I1 or n °t) and unstable at threshold (for (3 SS = e c /la)- in the 
special case that 71 = 72 = 73 = 7 we may find simple analytical solutions as 

Al,2 = 7o, 
^3,4 = 1, 



e 



A 5 ,6 = 7 + — yxi-xl 
7o 



A 7 ,s = 1--^X\-Xl (19) 

7o 

It is immediately obvious that only X 78 can possibly have a negative real part in this special 
case, which will happen when 

e 2 > -^L (20) 
X1-X2 

and means that any fluctuations will tend to grow, invalidating any linearised fluctuation 
analysis in this regime. We see that this is consistent with the critical pump value given 
above, in Eq. [131 We find numerically that the fluctuations above threshold cannot be 
linearised due to the presence of a zero eigenvalue. 

In Fig. CD, we give a plot of the different stability regions, for 70 = 71 = 73 = 1, 72 = 37i, 
and Xi = O.OI71, as xi an d e are varied. We see that when the system is oscillating and the 
non-pump modes are occupied, it is always unstable and must be treated numerically using 
stochastic equations. In the stable region below and to the right of the solid line we may 
use a linearised fluctuation analysis to calculate the correlations of interest. 

IV. DETECTION OF TRIPARTITE ENTANGLEMENT 

There are a number of inequalities whose violation is sufficient to demonstrate the exis- 
tence of continuous- variable tripartite entanglement, all of which are based on the insepa- 
rability of the system density matrix. Unlike bipartite entanglement, where two modes are 
either entangled or not, there are a number of cases to be considered, depending on possible 



partitions of the density matrix 30j. In this work we are interested in the case where the 
density matrix is not separable in any form, often known as genuine tripartite entanglement. 
Before we describe the criteria we will use here, we need to define the quadrature operators 
we will use, as different normalisations exist in the literature and can alter the exact form 
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FIG. 1: (colour online) Stability of the steady state solutions with xi = 0.01, 7o = 7i = 73 = 1, 
and 72 = 3, as \2 and the pump amplitude are varied. The dashed line shows the separation 
between the system with and without threshold. 

of the inequalities used. As we are considering that the cavity will be at resonance for all 
modes, we may use the orthogonal quadrature definitions, 

Xj = (Lj + Cbj, 

y. = -j^.-aj), (21) 

with the Heisenberg uncertainty principle requiring that V(Xi)V(Yi) > 1. We note here 
that any cavity detuning or Kerr interaction can change the quadrature angle at which the 
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321 ] , but this is not generally the case for a resonant 



best quantum correlations are found 
cavity with interactions. 

We will use two different sets of conditions to investigate the presence of entanglement 



in this system, both of which were described by van Loock and Furusawa 
these gives a set of inequalities, 

V 12 = ViX, - X 2 ) + V{Y X + % + g 3 Y 3 ) > 4, 

10 



The first of 



V 13 = ViX, - X 3 ) + V{Y t + g 2 Y 2 + Y 3 ) > 4, 

V 23 = V(X 2 - X 3 ) + V{g x Y x + % + Y 3 ) > 4, (22) 

the violation of any two of which shows that the system is fully inseparable and genuine 
tripartite entanglement in guaranteed. The are arbitrary real numbers which may be 
chosen to minimise the correlations, and will be optimised here as was done in Ref. 
giving 

V{Y u %) + V{Y l: Y 3 ) 
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92 



V{Y 1 ) 
V{Y X) Y 2 ) + V{Y 2 ,%) 



2) 



= _v(y,.y,) + v(r 3 ,r,) 3 

V(Y 3 ) 

The second conditions provide inequalities for which, if any one is violated, genuine 
tripartite entanglement is demonstrated. They are 

Vl23 = vfr - *±^) + V(Y t + > 4, 

V 3U = V(X 3 - ^Jl) + V(% + ^) > 4, 

V 23 i = V(X 2 - + V(% + > 4. (24) 

As in previous cases where the system is described by an asymmetric Hamiltonian (i.e. mode 
indices cannot be swapped without changing the system), the correct choice of indices during 



the measurement is important for both sets of correlations given above [151 ]. 

All these correlations can be simply calculated from the intracavity spectral matrix of 
Eq. [11] and the use of the standard input-output relations [2l| to give the measurable spectra 
outside the cavity. For example, spectral variances and covariances are calculated as 

S^(u) = 1 + 2 1j S Xj (u), 

SxUM = Zy/WiSxiM"), (25) 

and similarly for the Y quadratures. As this notation can become rather clumsy, we will 
use Sij and Syk in what follows to refer to the output spectral qualities equivalent to the 
Vij and Vijk correlations defined above. The same inequalities hold for these. 
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V. SPECTRAL RESULTS IN THE STABLE REGIME 



Although it is possible to obtain analytical results for the SV,- and Sijk, these are extremely 
unwieldy and not at all enlightening. We have therefore chosen to present the results graph- 
ically for various parameter regimes. In Fig. [2] and Fig. [31 we show the results for the two 
different types of correlations at half the critical pumping amplitude, in the regime where 
X2 = 0.4%i and with the loss rates set as 70 = 71 = 73 = 1, and 7 2 = 3. In all results 
presented here we have used a value of xi — 0.01. In Fig. [2] we see clear evidence of genuine 
tripartite entanglement, with both S123 and S312 obviously violating the inequality, with only 
one of these being below 4 already being sufficient. In Fig. [3J where two of the inequalities 
need to be violated, we see that an d both show entanglement, although not over as 
large a frequency region as the Sijk- As in Ref. 1||, this is a result of the asymmetry of the 
Hamiltonian and the fact that a violation of the tripartite inequalities is a sufficient but not 
necessary condition for the demonstration of tripartite entanglement. 




co (units of y ) 



FIG. 2: (colour online) Sijk criteria below threshold for e = 0.5e c , 70 = 71 = 73 = 1, and 72 = 3. 
In this and all subsequent graphs, the results are dimensionless. 



We see in Fig. H] and Fig. that, with e = 0.9e c , the violation of the inequalities has 
increased for two of the S^k, with the spectra bifurcating so that no entanglement is seen at 
near zero frequency. The Sij also show increased violation as the threshold is approached, 
but again not near zero frequency. The again do not indicate full inseparability over as 
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FIG. 3: (colour online) S%j criteria below threshold for e = 0.5e c , 7o = 7i = 73 = 1, and 72 = 3. 

wide a frequency range as the Sijk- 

In the region without an oscillation threshold, that is X2 > vTipYiXii we m ay also apply 
the linearised analysis. In Fig.[6]and Fig. Owe show the correlation functions for \2 — 2.5xi, 
with e = 1.5e° and the other parameters unchanged from Fig. [2j We see that entanglement 
is found in this regime and that its demonstration is less dependent on which particular 
correlations are measured, although 5*312 and 5*23 do not violate the inequality by a large 
amount. The main conclusion to be drawn from these results is that the S^k correlations 
are the most appropriate to use for this system, with S123 giving the maximum violation of 
the inequalities. 

We will now investigate the effects of changing e on the correlations, concentrating on 
the Sijk, as these have proven to be a more sensitive measure than the Sy. We firstly 
examine the region which has a threshold, that is where e < e c and X2 < Xi \J 12/11- We 
will present these results at the frequency of maximum violation of the inequalities, with 
the range < uo < 70, rather than fixing the frequency as the other parameters are changed. 
In Fig. [8] we show the minima of the three Syk as the pump varies between zero and the 
critical value, with xi — 0.01, X2 = 0.4^1, 7o = 7i = 73 = 1, and 72 = 3. We see that 
in no case does S231 violate the inequality, while the other two show clear violations, with 
S123 decreasing as threshold is approached. We note here that the results in the immediate 
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FIG. 4: (colour online) criteria below threshold with e = 0.9e c 




co (units of y ) 

FIG. 5: (colour online) Sij criteria below threshold with e = 0.9e c . 

neighbourhood of the threshold are not expected to be accurate, due to the invalidity of the 
linearised fluctuation analysis at that point. 

In the parameter regime where there is no oscillation threshold, we can also investigate 
the effects of varying both the pumping and \2- We again set 7 = 71 = 73 = 1, 72 = 3 
and Xi = 0.01, and will allow uj to vary so as to find the maximal violations. In Fig. [9] and 
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FIG. 6: (colour online) Sijk criteria for the system without threshold, with e = 1.5e°. 




FIG. 7: (colour online) Sij criteria for the system without threshold, with e = 1.5e°. 

Fig. [10], we plot the S^k as a function of e/e°, for values of X2 = 2%i and 3xi- We again see 
that 5*123 gives the maximal violations, although this does not increase monotonically with 
pump amplitude. In Fig. [TTJ and Fig. [12] we show how the correlations, again at the optimal 
frequencies, change as \2 is increased from x^"* ( = Xi\H^/7i)> f° r ^ wo different pumping 
amplitudes. We again see clear evidence of genuine tripartite entanglement over the range 
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FIG. 8: (colour online) Minimum of the Sijk at any frequency between and IO70 as the pump 
varies up to e c in the region with threshold. 

shown, with S123 again showing the maximum violations of the inequality. We note that the 
violations do not increase as \2 increases, but that 6*123 has its minima at y^ lt . 




FIG. 9: (colour online) Tripartite entanglement criteria for the system without threshold, with 
70 = 71 = 73 = 1, 72 = 3, xi = 0-01, and X2 = 2*1 . 



The results found from our semi-classical analysis show a range of different behaviours, 
including threshold behaviour which depends on an interplay of the two nonlinearities and 
could not be found in the undepleted pump treatment of Ref. [19l |. We have been able to 
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FIG. 10: (colour online) Tripartite entanglement criteria for the system without threshold, with 
70 = 71 = 73 = 1, 72 = 3, xi = 0.01, and X2 = 3xi- 

demonstrate that genuine tripartite entanglement is available for a wide range of parame- 
ters but have not been able to analyse the system in the above threshold regime, due to 
phase diffusion of the modes and the inapplicability of the linearisation procedure. In or- 
der to investigate the behaviour above threshold we will now turn to numerical stochastic 
integration. 



VI. STOCHASTIC INTEGRATION IN THE UNSTABLE REGIME 



In this case the numerical integration of the full positive-P equations (Eq. [6]) presents 
stability problems, so that, although they were useful in deriving the correct linearised equa- 
tions to calculate normally-ordered correlation functions, we will now turn to what is known 
as the truncated Wigner representation [28[]. Following the standard procedures [231 ]. the 
generalised Fokker-Plank equation for the Wigner representation pseudo-probability func- 
tion of the system is found as 

dW { 
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FIG. 11: (colour online) Tripartite entanglement criteria for the system without threshold, with 
7o = 7i = 73 = 1) 72 = 3, Xl = 0.01, and e = 0.5e°. The results are plotted as a function of 
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FIG. 12: (colour online) Tripartite entanglement criteria for the system without threshold, with 
7o = 7i = 73 = 1) 72 = 3, xi = 0.01, and e = 0.9e°. The results are plotted as a function of 
Xf* = Xi a/72/7i- 
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We immediately see that the above equation contains third-order derivatives so that it 
cannot be mapped onto a set of stochastic differential equations. Although methods have 
been developed to map these type of generalised Fokker-Planck equations onto stochastic 
difference equations in a doubled phase space [34 1 , the integration of these can present more 
stability problems than the positive-P representation, so we will not pursue this approach 
here. Hence we neglect the third-order derivatives to allow a mapping onto the set of 
stochastic equations in the truncated Wigner representation, 

dai 



dt 
da 2 

~dt 

das 

d/3 



7iai + Xi^lP + d — (m + iV2), 

72«2 + X2«3/? + d — [V3 + m), 

/7i" 

73«3 + Xi<*iP ~ X2a 2 p* + d — {r] b + irje), 
X2a 2 a* 3 + \p^(m + im) 



Xl«l«3 



(27) 



where the r)j are Gaussian random noises as defined by Eq. [7J As well as not containing 
multiplicative noise terms, another important difference from the positive-P equations is 
that the initial conditions on each stochastic trajectory must be drawn from the appropriate 
Wigner distribution for the desired quantum state of the mode. We will be beginning our tra- 
jectories with vacuum inside the cavity, so that, for example, we choose a™(0) = (£f 
on the nth trajectory (and similarly for P(0)), where the £ are normal Gaussian random 
numbers with zero mean. The Wigner representation naturally calculates symmetrically- 
ordered operator averages, so that care must be taken with any necessary reordering to give 
predictions for observables. We also note here that, while cases have been found where 
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the truncated Wigner representation can give inaccurate results |35j, |36j, we expect it to be 
accurate here because all four modes are macroscopically occupied in the region that we 
are using it to investigate. In our stochastic integration we have set 7o = 7i = 73 = 1, 
72 = 37i, Xi = O.OI71, X2 — 0.4%i and e = 1.5e c . An indication of the accuracy is that it 
gives predictions for the intracavity field intensities that are consistent with the analytical 
values given above, as shown in the table. 





Analytic 


Wigner 


\(3\ 2 


1.056 x 10 4 


1.056 x 10 4 


cti 2 


5.0143 x 10 3 


5.0141 x 10 3 


\a 2 2 


89.1424 


89.0909 


l«a| 2 


4.7468 x 10 3 


4.7471 x 10 3 



In Fig. [13] and Fig. [THwe show the intensity of the a modes inside the cavity, demonstrat- 
ing both that we have reached the steady-state regime and that the modes are macroscop- 
ically occupied. Integration for values of e close to e c typically took much longer to reach 
the steady-state, due to critical slowing down, a well-known phenomenon which occurs in 
the vicinity of phase transitions. Due to the phase diffusion predicted in the analytical so- 
lutions, the averages of the atj are essentially zero and no entanglement is registered by the 
Vijk correlations in the steady-state, with these all being far from violating the inequalities. 
As shown in Fig. [151 there is some violation in the initial transient regime before the pump 
mode within the cavity builds up to its threshold value and the system begins to oscillate. 
This transient feature is unlikely to be of any practical use, as it exists for only a few cavity 
lifetimes and the fields are no more intense than in the below threshold regime, where gen- 
uine tripartite entanglement is readily seen in the steady-state regime. We note here that 
the three-mode EPR correlations [3] give much smaller values in the steady-state, but still 
do not violate the inequalities. It is possible that a small injected signal at uj 2 could serve to 
lock the phases and enable entanglement to be observed, but investigation of this is outside 
the scope of the present work. 

VII. CONCLUSION 



We have analysed the system of singly-pumped intracavity coupled downconversion and 
sum-frequency generation with a quantised pump field. Unlike previous analyses, this en- 
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FIG. 13: (colour online) Intensity of the output modes |«i| 2 , obtained by stochastic integration of 
the truncated Wigner equations averaged over 2.65 x 10 5 trajectories. Note that the horizontal 
axis is now time scaled in units of 7j~ • 

abled us to define the threshold properties of the system and analyse the dynamics with 
all modes oscillating macroscopically inside the cavity. One of the features we have found 
is that for some values of the experimental parameters, the threshold value of the pump 
field diverges so that, however strongly the cavity is pumped, the system will not oscillate. 
We have found that genuine tripartite entanglement is available in both the regions below 
and without threshold, but that, as in other systems with asymmetric Hamiltonians, not 
all measurable correlations will detect the violation of the entanglement inequalities. Above 
threshold the converted modes undergo phase diffusion, which prevents the detection of en- 
tanglement based on quadrature measurements except in the early transient regime. This 
signifies that the system is not a good candidate for the production of bright entangled 
output beams, unless a method can be found to overcome the problem of phase diffusion. 
However, it is still useful for the production of genuine polychromatic tripartite entanglement 
in all except the above threshold regime. 
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FIG. 14: (colour online) The Vijk for e = 1.5e c , obtained by stochastic integration of the truncated 
Wigner equations averaged over 2.65 x 10 5 trajectories. Note that the horizontal axis is now time 
scaled in units of 7-f 1 . 
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